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Abstract 

At tree- level Higgs production in association with a 6-quark pair proceeds through 
the small Yukawa bottom coupling in the Standard Model. Even in the limit where 
this coupling vanishes, electroweak one-loop effects, through the top-Higgs Yukawa 
coupling in particular, can still trigger this reaction. This contribution is small for 
Higgs masses around 120GeV but it quickly picks up for higher Higgs masses es- 
pecially because the one-loop amplitude develops a leading Landau singularity and 
new thresholds open up. These effects can be viewed as the production of a pair 
of top quarks which rescatter to give rise to Higgs production through WW fusion. 
We study the leading Landau singularity in detail. Since this singularity is not 
integrable when the one-loop amplitude is squared, we regulate the cross section by 
taking into account the width of the internal top and W particles. This requires 
that we extend the usual box one-loop function to the case of complex masses. We 
show how this can be implemented analytically in our case. We study in some detail 
the cross section at the LHC as a function of the Higgs mass and show how some 
distributions can be drastically affected compared to the tree- level result. 



1 Introduction 



The LHC will soon start running and collecting data. Although one expects surprises, 
discovering the Higgs is the highest priority. A lot of effort has gone in calculating the 
rate of production of this particle, within the Standard Model and beyond, for a host of 
channels and signatures, see [H [2] for a review. 

Higgs production in association with a pair of bottom quarks is not, especially in the 
Standard Model, a discovery channel since the coupling of the Higgs to the bottom quark 
is given by the small, Oim^jv) bottom-Higgs Yukawa coupling, where m;, is the bottom 
mass and v ~ 246GeV the scale of electroweak symmetry breaking. Nonetheless given the 
special role that can play the third generation of fermions in the mechanism of symmetry 
breaking and in particular the top-bottom quark doublet, a reconstruction of this Higgs 
coupling to bottom quarks is important. This reconstruction and interpretation of the 
measurements requires theoretical predictions that go beyond the tree-level approxima- 
tion. Many of these calculations, most of which concern the important QCD corrections, 
have already been performed[3]. Usually one expects the electroweak corrections to be 
small and not compete with the QCD corrections. However, one should bear in mind 
that the top Yukawa coupling Oim t jv) is of order the strong coupling constant. If this 
coupling takes part in the electroweak corrections the latter may not necessarily be small. 
Other Yukawa couplings that are not negligible are the Higgs Yukawa couplings]. Both 
these couplings are involved when one considers the electroweak corrections to bbH pro- 
duction at the LHC. Another important property of the electroweak effects is that this 
cross section can be triggered off by one-loop corrections involving the top quark and W 
gauge boson (or Goldstone) loops even for vanishing bbH (or m& = 0) coupling, where the 
Born cross section vanishes. 

We[l] have, very recently, studied the effects of the leading (Yukawa-type) electroweak 
corrections for bbH production at the LHC in a situation where both 6's are tagged, re- 
quiring somewhat large pt b, as would be relevant for a measurement of the bbH couplings 
and a complete identification of this channel. The study we performed concentrated on a 
Higgs with a mass below 150GeV not only because this range is preferred by the precision 
electroweak data but also because the cross section decreases much with increasing Higgs 
mass. It was found that, after all, the Next-to-Leading Order (NLO) corrections were 
small and could be safely neglected. In the limit where the bbH coupling vanishes and 
where the cross section is induced solely through electroweak loops, we found that this 
effect is much larger than the NLO correction and increases rapidly with the Higgs mass. 
We pointed out that, for this contribution, as Mh > 2Mw our perturbative calculation 
becomes unreliable since the loop integrals start showing numerical instabilities. We had 
identified this behaviour as a leading Landau singularity (LLS)[SJ[n] which is a pinch sin- 
gularity of the loop integral. This, in part, has an interesting physical origin: the on-shell 
production and rescattering of the top quarks into on-shell W bosons, the latter giving 
rise to Higgs production through WW boson fusion. This LLS of the one-loop four-point 
function is not integrable when one considers the square of the loop amplitude as needed 

*As this paper is on the Yukawa corrections neglecting corrections of order the electroweak gauge 
coupling , we use the terminology Higgs Yukawa coupling for the Higgs self-coupling which in the Standard 
Model is not a gauge coupling. 
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for vanishing bbH coupling. The NLO contribution, on the other hand, is integrable. 

The aim of this paper is to extend the study we made in [I] to higher Higgs masses. 
The emphasis will be on the LLS problem and the pure one-loop contribution in the 
limit of vanishing m& since this is the major hurdle. For completeness we will also give 
results for the NLO contribution for Higgs masses not covered in our previous calcula- 
tion. Beyond the phenomenological impact of the LLS for the case at hand, the study 
of the LLS in this paper should be of interest for other situations considering that one 
rarely encounters such singularities, as compared to the inverse (vanishing) Gram deter- 
minant which is not a genuine physical singularity but an artifact of the reduction of the 
tensorial integrals. Some of the few examples in the relatively recent literature where 
some aspect of a Landau singularity shows up include ZZ — > ZZ^7\ and the 6-photon 
amplitude |H] in the Standard Model both with massless particles in the internal states 
involving a four-point function. Beyond the Standard Model we can mention loop cor- 
rections to sfermion pair production in supersymmetry[9| and Higgs production from the 
decay of a fourth generation b-like quark [in], both these examples involve heavy instable 
particles in a three-point function. In [9] no special treatment of the singularity is required 
since the study is made at the NLO level where this singularity is integrable. In[T0] the 
width of the internal unstable particle is called for. In ZZ — > ZZ, the study[7] keeps 
away from the region of the LLS, while it is argued that the LLS should disappear if one 
considers a more inclusive cross section where the Z boson would decay or the initial Z 
are grafted to light stable fermions. For the case of the 6-photon amplitude the situation 
is quite subtle. The QED dynamics is such that the LLS disappears at the level of the 
total gauge invariant amplitude after summing on individual diagrams. The LLS issue 
can also be relevant for the nascent cut techniques of computing loop amplitudes, for a 
recent review seepj]. This is the reason we devote a fair part of this paper to the study 
and solution of the LLS. Our solution to the problem of the LLS for Higgs production 
through gg — > bbH is to endow the resonating internal particles, namely the top quark 
and W gauge boson with a width. The extension of the usual loop libraries, such as 
FF[T2"] of LoopTools[T2], to the case of complex masses is not trivial especially if one 
insists on an analytical implementation. We will show how the case at hand lends itself to 
a fairly manageable implementation of complex masses for the four-point function which 
is computer-time effective. The introduction of the widths avoids all numerical insta- 
bilities and smooths out the cross section when we enter the phase space region of the 
LLS. It rests that this effect can still give large corrections particularly for some specific 
distributions, like for example the pt distribution of the bottom quark or the Higgs boson. 

The plan of the paper is as follows. In the next section we set the framework for our 
calculation with a reminder on the SU(3) (QCD) gauge invariant classes of the electroweak 
contributions and the helicity properties of the amplitudes. We then briefly uncover the 
class and type of diagrams that contain a potential leading Landau singularity. Section 3 
follows with a general discussion on the Landau singularities first exposing the conditions 
under which such singularities can show up for the scalar iV-point function. We then 
carefully extract the nature of the singularity before moving into a detailed investigation 
of the scalar 4-point function at the origin of the LLS in our case, for gg — > bbH. Section 4 
discusses how this singularity can be regulated through taking into account the width of 



2 



the unstable particles running in the loop. Section 5 describes how these widths are im- 
plemented through a modification of the loop integrals that should be defined for complex 
masses of the loop particles. In particular we describe our analytical implementation of 
the complex masses suitable for our problem. We will also discuss the various checks we 
made to insure the correctness of the implementation. Section 6 gives briefly our input 
parameters and cuts and describe how the cross section at the pp level is obtained. Sec- 
tion 7 gives our main results for the cross section pp — > bbH at the LHC in the limit of 
vanishing Higgs coupling to 6-quarks. In this case the cross section is induced at one- 
loop and we need, in particular, to integrate the square of the 4-point loop integral over 
the kinematical phase space. This calls for our new implementation of the box one-loop 
functions. We will discuss the behaviour of the cross section as a function of the Higgs 
mass and study a few distributions. Section 8 turns to the NLO result for Mh > 150GeV, 
completing therefore the study we made in [I]. Section 9 summarises our findings. The 
paper contains also three appendices. In the first we give the details of our derivation of 
the nature of the singularity while the second appendix gives technical details about the 
handling of complex masses in one-loop scalar box functions. The third appendix details 
the singularities of the 3-point function. Many key issues about the LLS are unravelled 
in this case which help in better understanding the issues in the 4-point function. 

2 A quick reminder and general considerations of the 
one-loop electroweak structure 

At LHC energies the exclusive bbH production with both fe-quarks tagged is dominated, by 
far, by the gluon gluon initiated subprocess. We therefore only consider, as we have done in 
11], the gluon-gluon initiated subprocess g(pi, Xi)+g(p 2 , A 2 ) — > b(p 3 , A 3 )+6(p4, A 4 )+if(p 5 ). 
Aj = ± and Pi with i — 1, 2, 3, 4 stand for the helicity for the momentum of the particle. 
The corresponding helicity amplitude will be denoted as *A(Ai, A 2 ; A3, A4). 
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Figure 1: All the eight Feynman diagrams can be obtained by inserting the Higgs line to 
all possible positions in the bottom line. 



At tree-level the process is given by Higgs radiation off the 6-quark line, see Fig. [U 
The tree-level amplitude, -4o(Ai, A 2 ; A3, A4), is therefore proportional to XbbH the Higgs 
coupling to b. As has been done in previous analyses [T5J 021 H] , for the exclusive bbH 
final state, we will require the outgoing b and b to have transverse momenta |p^ 6 | > 20GeV 
and pseudo-rapidity \r] b,b \ < 2.5. These kinematical cuts reduce the total rate of the signal 
but also greatly reduce the QCD background. As pointed in [IT] these cuts also stabilise 
the scale dependence of the QCD NLO corrections compared to the case where no cut 
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is applied. In the approximation of neglecting the bottom mass the whole contribution 
vanishes, since the Higgs coupling to b vanishes. The massless bottom limit can also be 
taken, but by keeping XbbH as an independent parameter with a non zero value. In this 
limit the tree-level contribution consists of only the amplitude Ao(\\, ^2; A, — A)Q. This 
turns out to be a very good approximation with the cuts we have taken, seejl]. 

At the one-loop level the electroweak effects introduce a rich structure even in the limit 
where one takes the leading Yukawa (top and Higgs) couplings that are most easily given 
by the contribution of the top/charged Goldstones contribution in the Feynman gauge [I], 
see Fig. El At one-loop, the diagrams are classified into three QCD gauge invariant classes 
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Figure 2: All the diagrams in each group can be obtained by inserting the two gluon lines 
or one triple gluon vertex (shown in class (c)) to all possible positions in the generic 
bottom line, which is the first diagram on the left. 



as displayed in Fig. El The Higgs couples to the bottom quark in class (a), to the top 
quark in the class (b) and to the charged Goldstone boson in class (c). As shown in 
Fig. [2] each class can be efficiently reconstructed from the one- loop vertex bbH, depending 
on which leg one attaches the Higgs, by then grafting the gluons in all possible ways. 
The difference in the coupling structure is another indication that each group forms a 
QCD gauge independent subset, see jl] for details. The analysis of [1] reveals that the 
contribution of class (a) at NLO is about —0.1% and thus can be totally neglected. Class 
(a) contribution naturally vanishes in the limit A&6# = as does the tree-level. In this limit 
the process is loop induced and triggered by diagrams in classes (b) and (c). Moreover in 
the limit m b — > with A;,^ 7^ 0, the one- loop corrections induce new helicity structures 
compared to those found at tree- level in this limit. 

TThe helicity amplitude method and the convention we use in this paper for the definition of the 
helicity state are based on [H [13] . 
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When trying to extend the study we have performed in [1] for M# > 2My/ we encoun- 
tered severe numerical instabilities for the cross section involving the square of the one- 
loop induced amplitude, which is the only remaining contribution in the limit 0. 
At the level of the NLO which involves the interference term between the tree-level and 
one-loop amplitudes no instability was present. On close inspection it was found that 
the instabilities were only due to the contribution from class (c) in particular to the box 
diagrams, including the box obtained from the reduction of the pentagon diagrams as dis- 
played in Fig. H]). At the partonic gluon-gluon level it was found there is no instability for 
^/s~g~g~ < 2m t and that independently of Mh and s/^gg the result was completely stable for 
writ — Mw- These threshold conditions were a sign for the possible existence of a leading 
Landau singularity for the box diagrams whose square is not integrable. The pentagon 
diagram in class (c) has no LLS but contains a sub-leading Landau singularity which is 
exactly the same as the LLS of the box diagram, obtained through the reduction of the 
pentagon to boxes. Some triangle diagrams of class (c) have also LLS (see Appendix O) 
but they are integrable hence do not cause any numerical instability. Since such singu- 
larities are little known nowadays and hardly encountered though we have referred to a 
few examples from the relatively recent literature in the introduction, we will discuss the 
issue of the LLS, their location and the condition on their appearance in the next section. 

Before that, let us remind the reader that, to calculate the cross sections, we use 
the same helicity amplitude method as the one used and explained in [3]. Details of the 
renormalisation scheme, for the NLO, and the optimization implemented in our code are 
the same as in [3]. To check the amplitudes and the cross sections we perform (QCD) 
gauge invariance tests and verify that our results are ultraviolet finite, see [3] for details 
of implementing these checks. 

3 Landau singularities 

Part of the discussion in this section has been summarised in [TT] and relies on [HJ E] 
although a few results are new. 

3.1 Conditions for a Landau singularity and the nature of the 
singularity 

Consider the one-loop process -Fi(pi) +-^2(^2) + • • • F^(pn) — > 0, where Fj stands for either 
a scalar, fermion or vector field with momentum pi as in Fig. [3j The internal momentum 
for each propagator is qi with i — 1, . . . N. Each momentum g« is associated with one real 
Feynman parameter Xi respectively. The scalar N-point loop integral in D space-time 
dimension reads 








J (27r)^z£> 1 £> 2 ...iV 
q\ — m\ + is with e > 0, q,[ = q + and ^ = q*, 




(1) 
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Qi = Qt comes from the fact that the g-integration hypercontour is along the real axis, 
according to the (infinitesimal) is prescription. The Feynman parameter representation 
reads 

f°° N f d D a 1 

T- = T ( N)[ ^■■■■Mg-.-l) / {2w) % {xiDi + X2D2 + ... XnD .,^ P) 



i=l 



Because of the Dirac delta function, the integration boundary in the Feynman parameter 
space are Xj = 0, % — 1, . . . , N. Thus the only important condition on Xj is that they are 
real and not negative. The singularities are given by the Landau conditions [HJ [TSJ, [5] 

Vz Xi(q? -mf) = 0, 

Eii x i1i = °> ( 3 ) 
Qi = Qi- 

If Eq. ([3]) has a solution Xi > for every i G {1, . . . ,N}, i.e. all particles in the loop 
are simultaneously on-shell, then the integral has a leading Landau singularity (LLS). 
If a solution exists but with some x,- = while the other Xi S £1X6 positive, the Landau 
condition corresponds to a sub-leading Landau singularity. To keep the analysis general 
let us therefore assume that Eq. ([3]) admits a solution with X{ = for i — M + 1, . . . , N 
with 1 < M < N and Xi > for every i £ {1, . . . , M }. Eq. ([3]) would read 

Xi = for i = M + 1, . . .,N, 

qf = mf , Xi > for i — 1, . . . , M, (4) 

EM n 
i=1 x^i = 0. 

For M = N one has a leading singularity, otherwise if M < N this is a subleading 
singularity. Multiplying the third equation in Eq. (Tj0) by qj leads to a system of M 
equations 

QhXi + Qux 2 H QiaiXm = 0, 

Q2i£i + Q22X2 H Q2MXM =0, 

. (5) 

Qmi^i + <3m2^2 H QmmXm = 0, 
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where the Q matrix is defined as 



Qij = 2q i .q j = m i + rrij - (qi - qj) = m { + rrij - (n - rj) ; i, j 6 {1, 2, . . . , M}, (6) 

and use is made of the on-shell constraint,^, e. the second equation in (jl]). Note that in 
Eq. © Xi > 0. 

The necessary conditions for the appearance of a Landau singularity can be summa- 
rized as follows 

det(Q) = 

X 2 > 2 (7) 

qi = qt 

for i = 1, . . . , M. The last condition, already encoded in Eq. ([3]), will prove to be useful, 
as we shall see. 

It has been shown by Coleman and Norton [TH] that if the matrix has only one zero 
eigenvalue then these equations are necessary and sufficient conditions for the appearance 
of a singularity in the physical region. 

If some internal (external) particles are massless like in the case of six photon scattering [8], 
then some are zero, the above conditions can be easily checked. However, if the internal 
particles are massive then it is difficult to check the second condition in Eq. ([7]) explicitly, 
especially if M is large. In this case, we can rewrite the second condition as follows 

Xj = det(Q jM )/det(Q MM ) > 0, j = 1, . . . , M - 1, (8) 

where Qij is obtained from Q by discarding row i and column j from Q and det(QjM) = 
d[det(Q)]/(2dQ jM ), det(Q M M) = d[det(Q)]/dQ M M- If det(Q M M) = then condition 
Eq. (JED becomes det(Q iM ) = with j = 1, . . . , M - 1. 

The condition of vanishing Landau determinant means that Q has at least one zero 
eigenvalue. In general, Q has iV real eigenvalues Ai, . . . , A^. Consider the case where Q 
has only one (non degenerate) very small eigenvalue Aat 1, which is what is occurring 
in our present calculation for gg — > bbH. To leading order 

Av = — , o-i = AiA 2 . . . Aat„i 7^ 0, a = det(Q). (9) 

With V = {xi, x®, ■ ■ ■ , x° N } the eigenvector corresponding to A^, we define v 2 = V.V. 
We will assume that Aj > for i = 1, . . . , K and A^- < for j — K + 1, . . . , N — 1 with 
< K < N — 1. It can then be shown that in D-dimension (see Appendix 1X1) 



jdiv 



)Ar+i e ina Kv (4:7r) aD T(a D ) 



{T ° NU VT 2(^+3)/2 ^/(-l)2a Kai (l\ NV 2 - fc)c 

N-K+l N-D+l 

OiK = o a D = n • ( 10 ) 



This result holds provided a\ ^ or in other words that the matrix Q does not have a 
degenerate zero eigenvalue. A similar result for the nature of the singularity has been 
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derived in [18] in the general case of a multi-loop diagram including the behaviour of 
the non-leading singularity. The extraction of the overall, regular, factor which is the 
.fT-dependent part in Eq. (imp is more transparent in our derivation. As stressed earlier 
the above result holds provided a\ 7^ 0. This general result has been derived with the 
assumption that formally N— D+l > 0, however unlike in [18] we can trivially analytically 
continue the result by using dimensional regularisation with D = 4 — 2e so that we can 
easily derive the nature of the singularity from Eq. fllUp even for the case of iV < 3 in 
D = 4. For the box in 4-dimension, N = 4, D = 4, a — > and a± 7^ we get 

in(3-K)/2 

(To)div = ; (11) 

1 ° M V(-l) 3 -*det(Q 4 )-«r 1 ' 

This shows that (Tq is integrable but its square is not. In the case N = 3 (the triangle), 
D — 4, one gets (see Appendix |X] for an alternative derivation not based on dimensional 
regularisation but along the one followed in 



in(2-K)/2 

C^U« = - 77=^ m \ ln (^ 2 - ( 12 ) 

87r v /(-l) 2 - A A 1 A 2 

Tq and its square are therefore integrable. 

The situation becomes more complicated when Q has a degenerate zero eigenvalue 
which happens in the case of the box diagram obtained in the case of the 6 photon 
amplitude or gg —> W + W~ [20J with massless internal particles. In D = 4 and for N >6 
a leading Landau singularity does not obtain, see for example p. 115 of [6]. We leave 
some of these issues for another publication though and will concentrate here only on our 
process. 



3.2 Application to gg — ► bbH 

Having set the stage for the occurrence of the Landau singularities we now turn to check 
that the numerical instabilities found in gg — > bbH are indeed due a Landau singularity. 
We concentrate on the box diagram in Fig. H] which can contribute a leading Landau 
singularity. The leading singularity of the 3-point function relevant for our problem is 
studied in Appendix O and serves as good starting point for the discussion to follow. The 
associated 5— point function where both external gluons attach to the internal top quark 
has no leading Landau singularity but rather a sub-leading Landau singularity which is 
exactly the same as the leading singularity that appears in the box diagram in Fig. HI It 
is thus enough to study, in detail, the structure and the singularity behaviour of this box 
diagram. We will keep the bottom quark massless unless otherwise stated. 

Defining the invariants s = s gg = (pi +P2) 2 , s% = (p$ +P5) 2 , $2 = (p4 + P5) 2 , and the 
on-shell conditions pi = p\ = p\ = p\ = 0, pi = Mjj, the kinematically allowed phase 
space region leads to the constraint 

Ml<S!<s, Ml— < s 2 < Ml + s-si. (13) 

We need to keep these constraints in mind as the solution of the Landau equations 
may fall outside the phase space. 
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Figure 4: A box diagram contributing to gg — > 66if i/ioi can develop a Landau singularity 
for Mh > 2M^y and > 2m 4; ie. all the four particles in the loop can be simultaneously 
on-shell. 



In terms of these invariants, the scalar box integral depicted in Fig. H] writes, in the 
nomenclature of LoopTools for example, as 

T 4 ( Sl , s 2 ) = D (M 2 H , 0, s, 0, Slj s 2 , M$y, M*r, m 2 , m 2 ). (14) 

3.2.1 On-shell and real conditions on the internal momenta Oj 

For the leading Landau singularity of the box in Fig. HJ the on-shell conditions on the 
internal particles read as q\ = q\ = m\ = m\ = M^, ?| = q\ — tu\ = m 2 = m\. The 
condition of real ^ = (gf, qj means that q^ > 0. At each vertex, one has 

\{Ml ml m 2 i+1 ) = (m 2 - ( mi + m l+1 ) 2 ^j \M 2 - (m, - m l+1 ) 2 ^j > 0, 

M 2 = ( qi -q l+l ) 2 , (15) 

with the usual A kinematical function, vertex % is identified as the the vertex to which the 
vector qi points according to Fig. HI is the invariant mass of the external leg at vertex 
i. Applying the condition of Eq. (1151) for the cases i = 1,3 we get 



M H > 2M W and > 2m*. (16) 

This requires that the normal thresholds for top quark production and Higgs decay into 
a W pair be opened. 

Condition YliLi x i ( li = in Eq. ([3]) in the case of the leading Landau singularity with 
Xi > is nothing else but the addition of N vectors, XjOj, with norm |:r;mj|. This says 
for example that not all time components g° can be positive or negative. For N = 4, 
either one vector has sign(xjg?) opposite to all the other three or there are 2 vectors x^ 
with positive signs while the others have a negative sign. In our case it is easy to see 
that we can only take > 0,^2 3 < 0. These simple considerations furnish additional 
inequalities that are constraints on the appearance of a LLS. Applied at the four vertices, 
for example in the rest frame of one of the internal on-shell particle [21], these give the 
additional normal thresholds of this 4-point function 

m t > M w (17) 
si > {mt + M w ) 2 and s 2 > (m t + M w ) 2 . (18) 
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These strong requirements on the opening up of the normal thresholds will delimit the 
region where a LLS will occur, as given by the vanishing of the Landau determinant. 
These normal thresholds are also normal thresholds of the reduced diagrams, 3-point and 
2-point functions, obtained from Xi = and are necessary condition for a LLS for these 
integrals, see Appendix O 

The on-shell and real conditions on the internal momenta qi with Yl x i1i = 0' Xi > 
have been given a beautiful pictorial physical interpretation by Coleman and Norton [T9~|. 
Each qi can be regarded as the physical momentum of a physical particle, we can associate 
to the Feynman diagram a space-time graph of a process with on-shell classical particles 
moving forward in time, XiiUi can be regarded as the proper time of particle %. The 
vertices are regarded as space-time points. AXi = Xiqi (no sum over i) is a space-time 
separation. 



3.2.2 Landau determinant 



The necessary conditions given by the inequalities above having to do with the opening up 
of normal thresholds need to be supplemented by the requirements of a vanishing Landau 
determinant. The reduced matrix, which is equivalent in this case to the Q matrix 
for studying the Landau singularity, is given by 



/ 



S 4 



2A1, 



w 



2M w m t 
\ 2M w m t 



1 



2M w m t 



*2 



2M w m t 



2M w m t 

2Mw"mt 
1 



2M w m t 
™-t+ M w-s 2 
2M w m t 
2m? — s 



\ 



2m? 



2m\ 
1 



s l i 



Q 



2mirrij 



(19) 



/ 



With s and Mh fixed one can study the behaviour of the determinant as a function of the 
invariant s\ and S2- The determinant is a polynomial of order 2 in each of these variables. 
In terms of s% for example it reads 



det(Q 4 



16M^ v m t det(5 4 ) = as\ + bs 2 



'2; 



b 



A(s 1 ,m t 2 ,M 1 2 



2{-s\(m 2 t + M } 



si - (m t 
) + si[(m 



Mi 



w 



M 2 w f 



{S2 

2 }[si - (m t - M w ) 2 }, 

(s-2m 2 t )(M 2 H -2M 2 w )} 



-sM 



H 



m t + 
-M 2 W 



Ml)}, s° 2 = -b/2a, 
I 2 + 2M 2 H s{m 2 + M 2 v )s 1 



+ 



sM 2 H [{s - Am 2 t ){M 2 H - AM 2 W ) - 4(m 2 + M 2 W ) 2 } 
b 2 - 4ac 



4a 2 



(20) 



Writing det Q4 as perfect square in s 2) like above for example, and a remainder which 
is the discriminant of the quadratic form that does not depend on s 2 can be revealing. In 
our case we find 



det(Qzi 



- det Q 2 (si] m 2 , My/) [ (s 2 - s 



4) 2 



det g 3 (si, M 2 H , 0; m 2 , M w , M w ) det Q 3 {s u s, 0; M w , 



m 2 , m 2 ) 



det Q 2 ( Sl ;mj,M w ) 



delQ 2 ( Sl ;mlM w ) 



1(21) 
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det Q^s are the Landau determinants of the 3-point function sub-diagrams obtained from 
the original 4-point function by shrinking one internal line to a point, forming sub- 
diagrams where the invariant Si is an argument of these 3-point functions. Likewise 
for det Q 2 obtained by further shrinking one of the triangles. The corresponding 2- and 3- 
point functions are shown in Fig.0 Our convention for det (^3,2 as concerns its arguments 
is given in Appendix O The factorisation in Eq. (121j) can be derived [22] for symmetric 
matrices based on the Jacobi ratio theorem for determinants [23j. Each sub-determinant 
of the reduced three-point function can be further reduced into exactly such a factorised 
form, see Appendix O This makes the identification of the sub-leading singularities very 




Figure 5: Three reduced diagrams of the box diagram in Fig. [^|, that contain s\ as an 
invariant of the 2- and 3-point functions and whose Landau determinants are given in 
Eq. ( dip . The self-energy diagram has a normal threshold. The two triangle diagrams 
contain anomalous thresholds. Note that the singularity structure of the second diagram, 
in the s\,Mh variables, is the same as the triangle studied in much detail in Appendix 3 
but with Si — > s 2 , see Fig. [73 The singularities of the second triangle can be obtained 
from the first one by s ^ Mfj, m t <-> My/. 

transparent. For example, det Q 2 {s\] m t, My/) = corresponds to a normal threshold, see 
Eq. ( ITBl) . It occurs for y/sl = m t + Mw (a/sT = m t — Mw is outside the physical region 
for Higgs masses of interest). Obviously we could have written the quadratic form in any 
of the variables si, s 2 , M^, s, the completion of the determinant will be the product of the 
determinant of two sub-diagrams. 

3.2.3 Numerical investigation of the four-point function and the Landau de- 
terminant 

We will always take m t = 174 GeV and Mw = 80.3766 GeV. Our investigation starts 
by taking y/s = 353 GeV, Mh = 165 GeV. The behaviour of the Landau determinant, 
the real and imaginary parts of the 4— point function Tq are displayed in Fig. i as a 
function of Sx, s 2 within the phase space. We clearly see that the Landau determinant 
vanishes inside the phase space and leads to regions of severe instability in both the real 
and imaginary parts of the scalar integral. 

To investigate the structure of the singularities in more detail let us fix y^s~[ = 
\j2(ml + M w ) ~ 271.06 GeV, such that the properties are studied for the single variable 
s 2 . This will also exhibit the sub- leading Landau singularities related to the reduced di- 
agrams. In the variables s 2 these are exactly the same as the ones we uncovered through 
Eq. (12~T1) . They are represented in Fig. [Sallowing for si — > s 2 (and x 2 — > x±, £4 — > £3). 
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Figure 6: The Landau determinant as a function of s\ and s 2 (upper figure). The real 
and imaginary parts of D as a function of Si and s 2 . The figure for the real part of D 
has been rotated since the structure is best seen with this view. 



Fig. [7] is very educative. We see that there are four discontinuities in the function 
representing the real part of the scalar integral in the variable y/s2. 

• As s 2 increases we first encounter a discontinuity at the normal threshold J~s~2 — 



a/ = m t + Mw = 254.38 GeV, representing Hb —>■ Wt. This corresponds to the 
solution (for the Feynman parameters) 2^3 = and £2,4 > of the Landau equations 
and can be associated to a leading Landau singularity for the 2-point scalar integral. 

The second discontinuity occurs ^/s~2 = 257.09 GeV. This corresponds to an anoma- 
lous threshold of a reduced triangle diagram. This corresponds to the solution x 3 = 
and #1,2,4 > of the Landau equations, (see Fig. [5]). The singularity structure of 
this diagram is studied in more detail in Appendix O We can explicitly check that 
y/s~2 = 257.09 GeV corresponds to the condition of vanishing determinant. One solu- 
tion of this equation does not satisfy the sign condition, Eq. (jHJ) and is not even inside 
of phase space. As shown in Appendix O only one solution, see also Eq. fl0.78D . is 
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Figure 7: The imaginary, real parts of D and the Landau determinant as functions of 



acceptable with 

*f = ^[MUM 2 w + m 2 t )-{m 2 t -M 2 w )M H \lMl-AM^) (22) 




2(mf + ML) + (Ml - AM') 1 + 1 " w 1 - = 

^1-AML/Ml 

which gives y/s^ = 257.09 GeV. Note that one of the necessary conditions for this 
anomalous threshold to occur in the physical region is M# > 2M^. At this normal 
threshold the value of s% is st? = 2(m 2 t + ML), see Eq. (122)1 . 

The third discontinuity at y/s~2 = 259.58 GeV corresponds to the anomalous thresh- 
old of the reduced three point function obtained from the box diagram by contracting 



to a point the x\ line, see the third diagram in Fig. 0, so that detQs(si, s, 0; ML, 
mf, m 2 ) = 0. Analogously ^J~S2 = 259.58 GeV is given by 

4 = (*K 2 + i^Sr) - vV* ~ 4r^ 2 (m? - M 2 ,)) , (23) 



2{m 2 + M 2 w ) 




^l-Am 2 /s, 
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• The last singular discontinuity is the leading Landau singularity. The condition 
det(S < 4) = for the box has two solutions which numerically correspond to yfs^ = 
263.88 GeV or yfs^ = 279.18 GeV. Both values are inside the phase space, see Fig. [7J 
However after inspection of the corresponding sign condition, only yfs?, = 263.88 GeV 
(with X\ ~ 0.53, %2 ~ 0.75, x% ~ 0.77) qualifies as a leading Landau singularity. 
y/s^ = 279.18 GeV has x x » -0.74, x 2 « -0.75, x 3 « 1.07 and is outside the 
physical region. 

The nature of the LLS in Fig. [7] can be extracted by using the general formula (II lj) . 
With the input parameters given above, the Landau matrix has only one positive eigen- 
value at the leading singular point, i.e. K — 1. The leading singularity behaves aM 

Df v = 1 (24) 

16M2 / m2 v /det( 1 S 4 ) - ie 

When approaching the singularity from the left, det(S f 4) > 0, the real part turns singular. 
When we cross the leading singularity from the right, det(5*4) < 0, the imaginary part of 
the singularity switches on, while the real part vanishes. In this example, both the real 
and imaginary parts are singular because det(S , 4 ) changes sign when the leading singular 
point is crossed. 



3.2.4 The leading Landau singularity region in the (M H , y/s) plane 

In practice we will have to integrate over the s\ and s 2 variables to obtain the total cross 
section at the partonic level. We will also have to integrate over s = s gg to arrive at 
the cross section at the pp level. Moreover, we would like to study the behaviour of the 
cross section by varying Mh- It is therefore important to quickly localise the range or 
region in the (y/s, Mh) plane where the leading Landau singularity occurs. This approach 
should in fact be followed in more general cases to check if one might encounter a po- 
tential problem prior to carrying the full phase space integration procedure with the full 
matrix elements. Necessary (but not sufficient) conditions on Mh and y/s to have a LLS 
correspond to the opening of normal thresholds as given in Eq. (1161) . These are easy to 
guess and are contained in the last two equations of Eq. ([7]). We have however to solve 
all of Eq. (J7J) together with the constraint that one is inside the phase space Eq. (TT3l . 
This, in general, is too complicated to be done analytically in a situation like ours with 4 
variables (Mh, s, s 1: s 2 ) and 2 parameters (Myy, m t ). However numerically the algorithm 
that goes through all the conditions is quite simple to implement. For instance one can 
start with the Landau determinant written as a quadratic form in s 2 by first computing 
the discriminant of the quadratic equation and check whether the latter is positive or 
negative, assuming the solutions are in the physical region. If the discriminant is positive 
one checks if the corresponding solution does not conflict with the positivity solution as 
implemented in Eq. (j81) . If this condition is satisfied then there is a LLS. In our case the 
result is shown in Fig. [SJ We conclude that the LLS occurs when 2M\y < Mh < 211GeV 

■l- The singularity of the 3-point function is logarithmic, see Eq. (jT3J) . Fig. [7] shows two 3-point singu- 
larities which look as if better behaved within LoopTools. 
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Figure 8: The region of 



and 2m t < y/s < 457GeV. The range of the LLS region depends on Mw and m t . If 
mt /My/ < 1 then the first two conditions in Eq. (J7J) can never be satisfied. In particular, 
if m t /M\v = 1 then the Landau determinant can vanish but the sign condition cannot be 
realised. When M H > 210GeV or s/s > 456GeV the Landau determinant det(Qi) can 
vanish inside the phase space but the sign condition Xj > cannot be fulfilled. 

The region of the leading Landau singularity in Fig. [8] is a surface of singularities in the 
plane of the kinematical variables \fs = y/s gg , Mh- This is bounded by three curves. It is 
important to stress again that the horizontal and vertical lines or boundaries correspond 
to the normal thresholds. These lines are also tangent to the upper curve delimiting the 
surface of LLS. We will get back to this property later. 

The algorithm we have just outlined is very easy to implement. The importance of the 
sign condition is crucial in determining the boundary of the leading Landau singularity 
region which occurs when — > 0. We will come back to this point shortly. Before doing 
so, it is worth coming back to the behaviour of D as a function of S2 like what we have 
shown in (Fig. [7]) and see how the location of the leading Landau singularity and the 
other discontinuities (related to other thresholds) move as Mh is varied. 

As in (Fig. CD we fix y/s = 353GeV and Jsl but with y/sl = 260GeV for M H = 
159, 165, 190GeV. All the curves will therefore show the two-point function discontinuity 
(normal threshold) at \J = 254.38GeV and three-point function discontinuity a/s| = 
259.58GeV, see Eq. (1231) . The other three-point function discontinuity at y^F an d the 
leading Landau singularity, if at all there, will of course move. The results are shown in 
Fig. [91 

• For Mh = 159GeV only the normal threshold at \J , and the a/5| discontinuity 
show up as expected since Mh < 2Mw- 
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Figure 9: The real part of D as a function of y^ for various values of Mh- For 
Mh = 2Mw we have taken s± = 2{m\ + M^). For the other cases, we take si = 260 GeV. 



For Mh = 165GeV (M# > 2Mw), the other three-point singularity shows at 
257.09GeV, together with the LLS at y 7 ^ « 283.5GeV. As 



■> 2 ~ Mu.oucv . ™ y/sl is increased the 
LLS moves to smaller values of S2, closer to the three-point function singularity as 
can be seen by comparing with Fig. [7] for the same value of Mh but higher value of 
y/si. This will be a common feature with the other cases with Mh > 2Mw, till the 
LLS disappears from the physical region. For yjs{ < 260GeV, no LLS develops. We 
have the ordering yf~s\ 



,tw 

'2 



< 
b 2 



< \ S% < 



>LLS 



For M H « 173GeV, V '4 = V$2 = 254.38GeV i.e the s% threshold coincides 
with the normal threshold. The LLS starts showing up at ~ 274GeV when 
^ii = 260GeV and moves to smaller values of S2 as ^/si increases. We have the 

< ^ 



ordering ^/s 



tw 



< 



y 7 S2 LS . The coincidence y 7 



signals 



the termination a leading singularity in the 3-point function, see Appendix As 
we increase Mh the LLS moves to smaller values of S2 and the discontinuity 
disappears from the physical region. 

For the special value Mh = 190.88GeV, the singularity has moved out of the phys- 
ical region but now the LLS coincides with the location of the s| three-point function 

For M H > 190.88GeV the 



< 



LLS 



singularity. We therefore have y 7 
LLS disappears from the physical region. 

Finally, we consider the special case of the threshold Mh = 2Mw where 
271.06GeV. One has to change s\ in the range defined by Eq. (Tl3l) with the condition 
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Si > (mt + M\y) 2 to make the LLS appear. It is easy to find out that the LLS only 
occurs when ^fs{ = ^J2(m1 + M^) = 271.06GeV and the LLS position coincides 
with the position of the three-point singularity yfs^ • We have the ordering < 

For future reference, it is worth noting that the LLS region opens up rather sharply 
when the normal thresholds open up and the bulk of the region is concentrated around 
these thresholds. Already for Mh > 200GeV the region squeezes into a very thin line. 



3.2.5 The leading Landau singularity region: analytical insight 

We will take two approaches. The first one is based on the observation that the boundary 
of the singularity region corresponds to a coincidence of a leading Landau singularity 
with a sub-leading singularity, this is the termination of the LLS [Ml E2J [21] . The second 
approach starts directly from the constraint or equation given by the vanishing of the 
Landau determinant. The extrema of this equation with respect to a particular choice 
of kinematical variables will define the termination of the LLS. Interpreting the equation 
as that defining a surface or a hypercurve, the extrema are tangents to the surface and 
are parallel to the corresponding coordinate variables. This will become clearer when we 
expose the derivation. 

i) A study of the LLS in the 3-point scalar integral relevant to our problem is quite 
simple since this function does, for fixed m t , Mw, involve a very small number of variables. 
Yet the study, see Appendix O, reveals some very general features. There is an LLS region, 
or curve, that is bounded by the normal threshold. This a manifestation of the fact that 
at the boundary, the leading singularity moves to the sub-leading singularity [Ml E2J [21] . 
This is also a phenomenon we observed in section 13.2.41 

Let us now analytically derive the surface shown in Fig. [HJ or rather the curve representing 
its boundaries in the (Mh, s) range. The lower bounds are just given by the normal 
thresholds of the two-point function so that Mh > 2Mw and y/s > 2m t , see Eq. f fT6l) . 
For each value of (Mh, s) there is a curve of LLS's defined by ^(si, s 2 , \Mh, s) which is 
constrained by the vanishing of det Q^Si, s 2 ) and subject to the sign conditions. For this 
discussion about the (M H ,s) range it is sufficient to only keep the (si,s 2 ) dependence 
of detQ4- As we scan over (Mh,s) we span a surface of LLS's. The key observation is 
that the curves terminate at a point corresponding to a sub-leading singularity, in this 
case a leading singularity of one of the 3-point function sub-diagrams which itself will 
terminate at the 2-point singularity, i.e. the normal threshold. For instance, writing 
det(Q 4 ) as a quadratic polynomial of s 2 as we did in Eq. (|2"T1) . there are 2 three-point 
sub-LLSs given by each det Q3 in Eq. (12 ip vanishing. The solutions of the latter are 
given, respectively, by Eq. (1221) and Eq. (l2"3l) . Let us take for definiteness the sub- leading 
singularity corresponding to s 2 in Eq. (122D . The argument works just as well with the 
other 3-point singularity s 2 . The coincidence constraint implies, for s 2 for example, a 
solution s 2 = s 2 = S2 and detQ 4 (si, s 2 ) = (with the proviso about the sign condition). 
Exactly the same argument can be put but now solving for the variable Si and exploiting 
the fact that our problem is symmetric in s\ <-> s 2 . The coincidence problem or the 
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constraint we are looking for translates into 

s 2 = s 2 and det Q 4 (si, s 2 ) = 0, 
si = and det Q4(si, S2) = 0, =>- 
detQ 4 (s 2 ,s 2 ) = and s 2 = s 2 (25) 

Only one solution to det Q4(s 2 , s 2 ) = passes the LLS sign conditions, with 



s 2 = 2(m 2 t + M 2 W ) - sj(s - 4m 2 )(M 2 H - AM 2 ,). (26) 
Equating Eq. ( 1261) with Eq. [221 we arrive at the equation of the termination curve 

\J(s - 4 ™?) = ( M ^K - - K + K)\J(m 2 h - uq^j . (27) 

Observe that this equation shows, in a very transparent way, that all thresholds: 

m t > M w , M H > 2M W , > 2m t 

need to be open simultaneously. We can invert Eq. (|27p to write the solution in terms of 
Mh- To arrive at the same result, it is more judicious however to go through exactly the 
same steps but choosing s| instead of s^. We derive 

y/(M% - AM 2 W ) = -L {yS{m 2 - M 2 W ) - (m 2 + M 2 w )^j{s - 4m?)) . (28) 

The maximum value of Mh (s/s) is obtained by setting y/s = 2m t (Mh = 2Myy), i.e. 
when the LLS, the two 3-point sub-LLSs and the normal threshold coincide. We have 



AM 2 W <M 2 H < AM 2 W + (m ' M " w) 2 



Am 2 < s < 4m 2 + ^ ~ f w) . (29) 



(m 2 - M 2 " 
Ml 



or numerically, 

348.00GeV < yfs < 457.05GeV and 160.75GeV < M H < 211.13GeV. (30) 

Of course, these analytical formulae reproduce exactly the curve in Fig. [8] that was ob- 
tained numerically. For example, we have arrived at the same, unique, solution by taking 
s i,2 = $2 an d si ]2 = s| in turn. This also means that the curve is also given by 

s s 2 = . (31) 

This constraint gives directly the equation for the bounding curve and avoids having to 
solve for si or s 2 as is done as an intermediate step in Eq. f)26|) . 

ii) Another interesting interpretation of the bounding curve which also leads to Eq. (1311) 
is based on the following. The leading Landau singularity in the (s\, S2) plane is a solu- 
tion of detQ^si, S2) = supplemented by the sign conditions. With fixed values of the 
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internal masses, the constraint det Q^si, s 2 , s, M H ) = is a constraint on the kinematical 
invariants for which a LLS can occur. This therefore defines a surface of LLS singularities, 
which one may want to visualise in the plane (si, s 2 ) or (s, M H ). Within the plane (si, s 2 ), 
the extrema of this surface are given by the tangents to this surface which are parallel to 
the coordinate variables, in this case s\, s 2 [22], therefore 

b — = with detQ 4 (si, S2) = and 



ds 2 

d det Q 4 (si, s 2 ) 
dsi 



with detQ 4 (5i,s 2 ) = 0. (32) 



These conditions are best exploited by using the quadratic form of det Q^Si, s 2 ) in s 2 
(and Si) given in Eq. fl2"Tj) . The first equation in Eq. f[3"2"j) with the help of Eq. fl2T|) leads 
to 

detg 3 (si,M|)detQ 3 (si,s) = 0. (33) 

The second equation, using again the same quadratic form in Eq. ( 12TI) leads to 

d det Q 3 (s u Mjj) <9det Q 3 (si, s) 2 

det Q 3 (si, s) H q det Q 3 (si, M H ) = 0. (34) 

asi osi 

The constraints of Eqs. f)33|34l) then require either i) both sub-determinants in Eq. ( 12T1) 
to vanish, detQ 3 (si, M H ) = det Q 3 (si, s) = 0. The latter requirement is exactly the 
condition given in Eq. ([2"Tj) . The other solutions of Eqs. f)33ll34p give the boundaries 
related to the normal thresholds , ii) det Q 3 (si, M H ) = 9dct Q3(si,M H ) _ q implies 
see Eq. (1C.81I) the normal threshold Mh = 2Mw is reached , while the third solution 
iii) det Q 3 (.Si, s) = 9dctl ^^ Sl ^ = q corresponds to the normal threshold s = (2m t ) 2 . 
These equations for the boundary define the LLS region presented in Fig. [HJ Note that 
ii) and iii) can also be derived from i) if one insists on finding the extrema of the curve 
det Q 3 (si, Mfj) = for example. This is the same argument that is used in Appendix O for 
the three-point function. Here we can carry this argument one step further starting from 
the fact that det Q 3 = is a condition for the Landau singularity of a 3-point function. 
The extrema and tangent argument applied at this level will show that the range in 
Mh and s are given by the vanishing of the corresponding det Q 2 which give the normal 
thresholds, Mu = 2My/ and s = (2m t ) 2 . This derivation shows that when the normal 
threshold is met all singularities of the 2-, 3- and 4-point function coalesce. Observe that 
in Fig. [8] the lines given by Mh = 2M\y and = 2m t are not only boundaries of the 
LLS region but also tangents to the extremum bounding curve given by Eq. ( 1271) . 

The arguments given above can be applied to derive the bounding curve and the range 
of the LLS's in the (s\, s 2 ) plane after elimination of the variables {M H , s) and taking into 
account the normal threshold condition, si. 2 > (m t + My/) 2 as the lower bound. The 
starting point in this case is to express det as a quadratic polynomial in M H for 
example. The solution of the bounding curve is given by 

si~S2 = ^ 2 ~ ^| (^XisumlM^) - ^/\(s 2 ,mlM^ . (35) 
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This translates into the bounds 



{m t + M w f < si >2 < {m t + M w f + 



m t M w 

254.38GeV < JsH < 324.44GeV. 



numerically 



(36) 



4 The width as a regulator of the Landau singularity 

As we have seen the leading Landau singularity requires all internal particles to be on 
their mass shell, see for example Eq. (JTj). This is akin to the usual singularity that occurs 
on resonance for a massive particle. These equations also show that if any parameter rrii 
is complex with a non zero imaginary part, the singularity is avoided. For an unstable 
particle the width provides this imaginary part. As can be inferred from Eq. (|7j), math- 
ematically, the width effect is to move the Landau singularities into the complex plane, 
so they do not occur in the physical region (the real axis). For our problem, the Landau 
condition in the interpretation of Coleman and Norton through Eq. [TTJ, m t > Mw, clearly 
shows that the singularity develops because of the instability of the top quark. Therefore, 
in principle, one should only include the width of the top as a regulator. Including the 
width of an unstable particle, whereby the mass of the internal particle becomes complex 
effectively sums a subset of higher order Feynman diagrams thereby taming the Landau 
singularity [in]. On the other hand, if one goes to higher order to implement the width 
then we would not only induce a width for the top but also for the W. Therefore to be 
realistic one should include the widths of both the top quark, T t , as well as the width of 
the W, T w . 
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Figure 10: Effect of the width of the W , Ty/ and of the top, T t , on the real and imaginary 
part of the four-point scalar function. 

We take the simple prescription of a fixed width and make the substitution 

m\ -> m\ - im t T t , -> - iM w T w . (37) 

Applied to the case of our four-point function one sees in Fig. [10] that indeed the width 
regulates the LLS and gives a smooth result that nicely interpolates with the result at 
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zero width away from the singularity. The normal threshold and the 3-point sub-leading 
singularity are also softened. The real part of the 4-point function still shows a smooth 
valley at the location of the LLS after regularisation. For the imaginary part we note that 
after introducing the width the LLS singularity is drastically reduced with a contribution 
of the order of the sub-leading singularity. 

As we will explain in the next section and in more detail in Appendix B the introduction of 
the width in a four-point function requires careful extension of the usual 4-point function 
libraries. In the case at hand, as will be shown, the four-point function with complex 
internal masses can be written in an analytical form, albeit with a larger number of 
Spence functions compared to the case of real masses. 

In our calculation of Yukawa corrections where all the relevant couplings depend only 
on the top-quark mass, the Higgs mass and the vacuum expectation value v, we will keep 
m t , Mh and v real while applying rules (13 7p to all the loop integrals. 

One might ask whether the same prescription as in Eq. (1371) for the Higgs mass can be 
of any relevance. A justification for this will require to consider the corresponding process 
including the Higgs decays with among other contributions, "resonant contributions" with 
an integration over the propagator of the Higgs. At least on a diagram by diagram basis 
this will not solve the problem since for example one still has to deal with the same 4-point 
function but with Mfj replaced by a certain p 2 H , taking into account the fact the leading 
Landau singularity occurs for a wide range of Higgs masses and values of the invariant 
p 2 H . On the practical side, recall that compared to the top and W width of about 2GeV, 
for Higgs masses of about 2Myy the width of the Higgs is O.lGeV, more than an order of 
magnitude smallei[§. 

5 Implementation of complex masses in the loop in- 
tegrals 

We have implemented complex masses in all the loop integrals we encounter in calculating 
the cross section in the limit XbbH — where the tree-level prediction vanishes. In this 
limit we can also set the mass of the bottom quark to zero. In the S'?7(3)-gauge invariant 
classification of Fig [2J class (a) vanishes in this approximation. In fact we had shown 
[I] that even with mj = 4.62GeV class (a) is totally negligible. Although it is only class 
(c) that shows severe numerical instabilities due to the presence of a leading Landau 
singularity in the 4-point box function and non-leading singularity in the 5-point function 
we introduce the width in all diagrams of both classes (b) and (c). 

For the tensorial and scalar loop integrals with up-to three legs we rely on LoopTools 
[T3"] which handles complex masses in up to 3-point functions. The 5-point functions 
are reduced to 4-point functions according to [2SI EZj- The tensorial 4-point functions 

§Our calculation of the leading Yukawa effects involves the charged Goldstone boson in the Feynman 
gauge through which the W mass enters. One may question whether it is appropriate to introduce a width 
here for a Goldstone boson considering that a Goldstone is defined as a massless state. Independently of 
the width one should first question why the Goldstone has a mass here. The point is in any other gauge 
than the Feynman gauge we would have had to consider the effect of the Goldstone and W exchange to 
derive the leading Yukawa effects. The physical thresholds are therefore captured in the Feynman gauge. 
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are reduced to the scalar 4-point function and 3-point functions. We therefore have to 
calculate only the scalar 4-point function with complex masses. The analytical calculation 
of 4-point function with complex masses in the most general case is practically intractable. 
If one of the external particles is lightlike, the standard technique of 't Hooft and Veltman 
[28] brings some light although the result writes in terms of 72 Spence functions. In our 
example, gg — > bbH with massless bottom quarks, there are at least 2 lightlike external 
momenta in all boxes, including the ones derived from the pentagon diagrams. If the 
positions, in the box, of two lightlike momenta are opposite then we can write the result in 
terms of 32 Spence functions. If the two lightlike momenta are adjacent, the result contains 
60 Spence functions. The detailed derivation and results are given in Appendix [B] We 
have implemented those analytical formulae for the case of two massless external momenta 
into a code and added this into LoopTools 0. 

We have performed a variety of checks on the new loop integrals with complex inter- 
nal masses. First of all, for all the tensorial and scalar loop integrals (4- and 5- point 
functions), we have performed a trivial numerical consistency check making sure that as 
the numerical value of the widths is negligibly small, widths — > + , one recovers the well 
tested result with real internal masses. For the scalar loop integrals, the results are com- 
pared to the ones calculated numerically in the limit of large widths, e.g. T tt w — lOOGeV, 
we find an excellent agreement. Furthermore, for the scalar box integrals the results can 
be checked by using the segmentation technique described in [29|. The idea is the follow- 
ing. At the boundary of phase space where the Gram determinant vanishes, the 4— point 
function can be written as a sum of four 3-point functions. The 3-point functions with 
complex masses can be calculated by using LoopTools. In this way, we have verified 
with excellent precision that the results of the scalar 4-point functions are correct at the 
boundary of phase space. We have also carried out a comparison with a dedicated purely 
numerical approach based on an extension of the extrapolation technique |30j. We have 
found perfect agreement!!. 

In a second stage we have performed checks at the amplitude level. A very trivial 
one was to check that the results with the new loop library exactly match the ones with 
the standard loop library with real masses in the limit widths — > + . Another important 
check was to verify that the results calculated with complex internal masses are QCD 
gauge invariant, see [I] for this check. 

Since the leading Landau singularity is integrable at interference level, the NLO calcu- 
lation with X bbH 7^ performed in [I] can be trivially extended to the region of M H > 2M W 
by using the same method without introducing widths for unstable internal particles. 

1 /2 

However, there is a small problem related to the universal correction (SZjj — Sv) where 

1/2 

the wave function renormalisation of the Higgs oZjj related to the derivative of the Higgs 



two-point function becomes singular when Mh equal to 2Mw or IMv. |31j. We regularise 
this singularity by separately introducing the widths of the W and the Z. This singu- 
larity, contrary to the leading Landau singularity, is due to the Higgs being an external 
one-shell particle. Other ways for dealing with this problem have been discussed [32J. 



^The implementation for the case of one massless external momentum is straightforward. However, 
we have not done this yet since it is not necessary for our present calculation. 
II We thank F. Yuasa for sending us the results of the extrapolation technique. 



22 



6 Inputs parameters and kinematical cuts 

The input parameters are the same as given in [3]. We rewrite them here together with 
new inputs which are the widths of the unstable particles appearing in the calculation. 

a(0) = 1/137.03599911, a s (M z ) = 0.118, 

M w = 80.3766GeV(G / , = 1.16639 x 10~ 5 GeV" 2 ), M z = 91.1876GeV, 
m t = 174.0GeV, IV = 2.1GeV, Y z = 2.4952GeV, (38) 

the top-quark width is calculated at the tree level in the SM as 

87rV2mf 

where the bottom-quark mass has been neglected. The Cabibbo-Kobayashi-Maskawa 
parameter Vtb is set to be 1. Most of our discussion concerns the most interesting case 
of the limit XbbH ~> where as we have discussed at length, see also [3], the effect of the 
fe-quark mass other than in the Higgs coupling is totally negligible. Therefore we set this 
mass to zero when discussing this limit in section [71 For completeness we will also give 
results for the NLO corrections in section [HI which require XbbH 7^ 0. There we will set 
nib = 4.62GeV. When we refer to the leading order contribution we will have in mind 
the cross section at the Born level calculated with m& = 4.62GeV. The cross section from 
the one-loop amplitude squared with XbbH —> will, in a few instances, be normalised to 
this Born cross section to give a measure of the new electroweak effect and so as to allow 
comparison with the NLO corrections. 

We consider the case at the LHC where the pp center of mass energy is yfs = 14TeV. 
Neglecting the small light quark initiated contribution, see [3], we use CTEQ6L [33 l 1341, 1351 
[36] for the gluon density function in the proton. The factorisation scale for this density 
and the energy scale for the strong coupling constant are both chosen to be Q = M z for 
simplicity. 

As has been done in previous analyses [TS1 EH [3], for the exclusive bbH final state, 
we require the outgoing b and b to have high transverse momenta |p^ 6 | > 20GeV and 
pseudo-rapidity \rj b,b \ < 2.5. These kinematical cuts reduce the total rate of the signal 
but also greatly reduce the QCD background. As pointed in [17] these cuts also stabilise 
the scale dependence of the QCD NLO corrections compared to the case where no cut is 
applied. In the following, these kinematical cuts are always applied. 

7 Results in the limit of vanishing X^h 
7.1 Total cross section 

We start with the cross section in the case where X^h = 0. In [3] we reported on 
results up to M H = 150GeV that showed that this cross section was rising fast as one 
approached the threshold M H = 2M W . Beyond this threshold our integrated cross sections 
showed large instabilities. As we discussed in section [3] this is due to the appearance of 
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a leading singularity which as we have advocated can be cured by the introduction of 
a width for the unstable top quark and W gauge boson. We also showed in section [3] 
that the region of Landau singularity spans the region 2M\y < Mh < 211GeV with 
2m t < y/s^ = \fs < 457 GeV, see Fig. El Before convoluting with the gluon distribution 
let us briefly look at the behaviour of the partonic cross section gg — > bbH paying a 
particular attention to this leading Landau singularity region. 




M H [GeV] Vs[GeV] 

Figure 11: Left: the cross section for the subprocess gg — > bbH as functions of Mh for 
various values of y/s including the case \/s = 2m t = 348GeV . Right: the cross section 
for the subprocess gg —>■ bbH as functions of y/s for various values of Mh including the 
case M H = 2M W = 160.7532GW. 



Figs [TTI show that indeed the widths do regulate the cross section. Moreover it is within 
this range that the cross section is largest even after being regulated. The (highest) peak 
of the cross section occurs for a Higgs mass of 163GeV about above the Mh = 2Mw 
threshold and for y/s = 351GeV about 2T t above the y/s = 2m t threshold. Figs. [TTI 
show that the cross section exhibits a peak structure close to the onset of the normal 
thresholds in Mh,\/s even when one is slightly outside the leading Landau singularity 
region of the 4-point function. In fact, this enhancement at the normal threshold is far 
from being totally due the 4-point LLS especially after the latter has been regularised by 
the introduction of the width. At the normal threshold there is an enhancement from the 
accumulation of all the 2-point, 3-point and of course the 4-point function. Moreover as 
we noted in section HI see Fig. QJJ], the introduction of the widths brings the contribution 
of the LLS to the level of a sub-leading singularity. 

The cross section at the pp level for the 14TeV centre of mass energy at the LHC 
as a function of the Higgs mass is shown in Fig. [12] taking into account the width of 
the top quark and the W gauge boson. For comparison we also show the cross section 
without the width effect outside the leading Landau singularity range of Mh- The sharp 
rise above Mh > 150GeV is nicely tamed. On the other hand note that on leaving 
the leading Landau singularity region around Mh = 211GeV, the width effect is much 
smaller and the figures suggest that one could have "entered this region from the right" 
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M H [GeV] M H [GeV] 

Figure 12: Left: the one-loop induced cross section as a function of Mh in the limit of 
vanishing bottom-Higgs Yukawa coupling for two cases: with and without widths. Right: 
the percentage correction of the contribution with widths relative to the tree level cross 
section calculated with \ bbH ^ 0. 



without having recourse to introducing a width. Indeed our numerical integration routine 
over phase space with the default LoopTools library does not show as bad behaviour 
until we venture around values of 2Mjy < Mh < 200GeV. The reason for this can be 
understood by taking a glance at Fig. [SJ For 200GeV < Mh < 211GeV the singularity 
region is considerably shrunk to a line so that one is integrating over an almost zero 
measure. The effect of the widths outside the singularity region is to reduce the cross 
section for M H = 120GeV, 140GeV and 150GeV by respectively 15%, 24% and 33% while 
for Mh = 210GeV, 230GeV and 250GeV the reduction is comparatively more modest 
with respectively 15%, 5% and 2%. 

Normalised to the Born cross section the new contribution represents a mere 2.6% 
for Mh = 120GeV. It increases however to as much as 49% for Mh = 163GeV before 
stabilizing to about 10% for larger Higgs masses. 

7.2 Distributions 

Effects of the new purely one-loop contribution being as large as ~ 50%, compared to the 
Born cross section even after being regulated through the introduction of the widths, it is 
essential that one looks at different distributions to see if this new effect can be described 
as a simple i^-factor. The two examples we show for Mh = 150GeV (before the onset 
of the leading Landau singularity) and for the Mh = 163GeV where the effect on the 
total cross section are largest show that the corrections are not uniformly distributed 
for all distributions. Figures [13] for Mh = 150GeV show the effect of the width. While 
the relative difference is rather uniform, about 33%, on the Higgs pseudorapidity, r] H , 
distribution, the transverse momentum distributions of the Higgs, p^, and the bottom p h T 
are strongly affected in particular for values which in the absence of the width showed a 
peak structure. There is still some peak structure in the pr distributions but the width 
effect reduces this by as much as 50%, while in the tails it is about 10%. 
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Let us now turn to M# = 163GeV. The correction, normalised to the Born cross 
section, for the Higgs pseudorapidity distribution is about 60% around the center region. 
The corrections to the px distributions can be enormous in some regions of phase space, up 
to 200% for the Higgs and about 170% for the bottom quark case. These huge corrections 
to the distributions in some region of phase space are again due to the effect of Landau 
singularities. 

One may question whether these large corrections signal the breakdown of perturbation 
theory and whether one expects (even) higher order effects to be large. We do not think 
so. First of all the relative large corrections have to do with the fact that for vanishing 
XbbH the tree-level cross section vanishes. Second, higher order effects have been captured 
in the introduction of the width and there is no reason to suspect that the leading Landau 
singularity we have encountered is affected by higher order effects. 

8 Results at NLO with X bbH ^ 

The results of the electroweak corrections at NLO which represent the interference con- 
tribution between the Born and the one-loop amplitude are much less interesting and 
numerically quite small, a trend that we had found already when studying at some length 
the electroweak NLO for Mh < 150GeV[lj. Moreover although some one-loop diagrams 
contain a leading Landau singularity at the interference level this singularity as we have 
shown in section [3] is integrable, see Eq. (TTTT) . The NLO contribution, apart from the 
Higgs wave- function renormalisation effect, is numerically stable even if one does not im- 
plement widths of the internal particles. The purpose of this section is to briefly present 
the results for the NLO. We first show that the effect of introducing the width is very 
small then show the NLO result without the internal widths being implemented hence 
these results are genuinely NLO results. These results thus complement the study we 
made for M H < 150GeV[l]. 

As discussed in section [2] the NLO Yukawa corrections consist of 3 QCD gauge invari- 
ant classes, see Fig. [2 Class (a) gives a totally negligible correction below 0.1%. We will 
not discuss this contribution any further here. Moreover, the leading Landau singularity 
we have discussed only shows up in class (c). As a first step we therefore study the NLO 
correction due to class (c) and weigh the effect of implementing the width of the inter- 
nal particles. Class (b) does not develop a leading Landau singularity and therefore the 
widths effects will be marginal. 

1 /2 

Another correction with enhanced Yukawa coupling is the universal correction, (SZg — 

1 /2 

5v) where SZ^ , the Higgs wave-function renormalisation constant involving the deriva- 
tive of the two-point function Higgs self-energy. The latter is ill-defined when Mh is equal 
to 2M\y or 2Mz- Here the width of all unstable particles, W, Z, t, will be kept o 
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Figure 13: The pseudo-rapidity of the Higgs and transverse momentum distributions of 
the Higgs and the bottom for Mh = 150GeV arising from the purely one-loop contribution 
in the limit of vanishing LO (XbbH — 0) f or two cases: with and without widths. The 
relative percentage contribution da(XbbH = 0)/daio is also shown. 
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Figure 14: The pseudo-rapidity of the Higgs and transverse momentum distributions of 
the Higgs and the bottom for Mh = lQ3GeV arising from the purely one-loop contribution 
in the limit of vanishing LO (XbbH — 0/ Its relative percentage contribution da(\f,bH — 
0)/daLo is also shown. 
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8.1 Width effect at NLO 

Our implementation of the width in the four-point function has been done in the limit 
of massless external quarks. To be fully consistent in the calculation of the one-loop 
amplitude with widths using the modified 4-point function we switch off the bottom mass 
in the spinors and propagators but keep XbbH ^ as an independent parameter. Our 
results for the NLO contribution of class (c) is shown in Fig. [15j First of all as we can see 
the overall correction is quite small, even at the onset of the (integrable) leading Landau 
singularity, the correction to the Born term is below 3.5%. The existence of a dip at the 
expected location is noticeable. Width effect softens the dip behaviour somehow but the 
effect is not as dramatic as what we have seen in the previous section for the loop squared 
results. We find that if M H < 158GeV or M H > 165GeV then the width effect change 
the NLO result but not more than 5% and are therefore totally negligible especially if one 
takes into account the smallness of the NLO result itself. Therefore the full NLO results 
can be studied by safely neglecting the width effect in classes (b) and (c). 

8.2 NLO corrections with rrib ^ 

The results for the NLO corrections are shown in Fig. [Tn]as a function of the Higgs mass. 
We implement widths only in the two-point function wave function renormalisation of the 

1/2 

**Note that SZ^ does not diverge when Mh = 2mt and the top-quark width thus has a marginal 
effect on SZli 2 . 
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Figure 16: Left: The relative NLO EW corrections normalized to the tree-level cross 
section, (b) and (c) correspond to the two classes of diagrams displayed in Fig. [H (8Z^ 2 — 
Sv) is the universal correction contained in the renormalization of the bbH vertex. "Total" 
refers to the sum of those 3 corrections. SZ 1 ^ 2 is calculated by taking into account the 
widths ofW, Z and the top quark. Right: the structure of (c)- correction which is a sum 
of two independent helicity configurations. 



Higgs. The latter contributes an almost constant — 1% correction apart from oscillations 
in the range 2Mjy to 2Mz due to the dips at 2Mw and 2Mz where the Higgs wave 
function is not analytic at those value. The effect of the widths of the W and Z smooths 
the behaviour and the corrections is never larger than 3.5% in this range of Higgs masses. 
The contribution from class (b) where the Higgs couples to the internal top decreases 
very slowly as the Higgs mass increases from HOGeV to 250GeV, as expected there is 
no structure as would be the case if this contribution were sensitive to any threshold or 
singularity. Class (c) on the other hand does, as expected, reveal some structure around 
Mh = 2Mw where we see a fall in the relative correction. The correction is however, 
despite this fall, quite modest ranging from ~ —1% for Mh = 160GeV to —4% for 
Mh = 210GeV. When we studied the effect of the width of the internal particles on class 
(c) at NLO, we did so in the massless quark limit. In that limit the outgoing quarks have 
opposite helicity so that only the 5a 3 ,-a 4 helicity amplitude survives, A 3i4 are the quark 
helicities. In our case here when the quark mass is reinstated, the <5a 3 ,a 4 helicity amplitude 
switches on. Fig. [TBT right) shows that these two helicity amplitudes behave differently as 
a function of the Higgs mass. The effect of the 6-quark mass makes the dip in the <5a 3 .-a 4 
much softer that in the massless case displayed in Fig. [15j In the <5a 3 ,a 4 the fall of the 
correction around Mh = 2M\y is more apparent. This is another manifestation of how 
the dynamics can affect the structure of a singularity. 

Adding the effect of all the contributions at NLO the total correction changes from 
-4% for M H = HOGeV to -8% at M H ~ 2M Z stabilising to around -7% past this value 
up to M H = 250GeV. 
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9 Conclusions 



At tree-level Higgs production in association with a 6-quark pair at the LHC is dominated 
by gg — > bbH where the Higgs is radiated from the 6-quark with a strength proportional to 
the bottom-Higgs Yukawa coupling. Unfortunately in the Standard Model this coupling is 
extremely small and therefore this mechanism is not a Higgs discovery channel although 
once the Higgs has been found the study of the Higgs coupling to the 6-quark through 
this reaction could probe interesting phenomena having to do with the the mechanism 
of symmetry breaking and the role played by the third generation fermions. Electroweak 
one-loop effects are usually small compared to the QCD corrections, however processes 
involving the bottom quark, electroweak one-loop corrections involve the top quark whose 
Yukawa coupling is of the order the QCD strength. More interesting for bbH production 
is that even in the limit where the bottom-Higgs Yukawa coupling vanishes and there- 
fore the Born tree-level cross section vanishes, electroweak one-loop effects, through the 
top-Higgs Yukawa coupling in particular, can still trigger this reaction. We studied these 
effects in some detail in a previous publication [1] but presented results for Higgs masses 
below 2M\y. We remarked that for the one-loop contribution in the limit of vanishing 
bottom-Higgs Yukawa coupling, the cross section was growing as the Higgs mass increased 
and that numerical results started showing instabilities past Mh > 2Mw- The aim of 
this paper was to extend the study performed in [I] to the mass range where numerical 
instabilities occurred. The origin of the numerical instabilities is due to the fact that 
some one-loop contributions, contained in some box diagrams, develop a leading Landau 
singularity. We have here reviewed in some detail the problem of the occurrence of the 
leading Landau singularity and investigated in more details the conditions and dynamics 
as concerns bbH production. Since this singularity is not integrable when the one-loop 
amplitude is squared, we regulate the cross section by taking into account the width of 
the internal top and W particles. This requires that we extend the usual box one-loop 
function to the case of complex masses. We show how this can be implemented analyt- 
ically in our case. We study in some detail the cross section at the LHC as a function 
of the Higgs mass and show how some distributions can be drastically affected compared 
to the tree-level result. For completeness we have also extended our study of the NLO 
Yukawa electroweak corrections which represent the interference between the one-loop 
amplitude and the tree-level amplitude. At this level the Landau singularity is integrable 
and therefore does not require that one endows the internal particle with a width. The 
NLO correction is found to be small. 
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Appendices 



A Nature of the leading Landau singularity 

We give in this section more detail about our derivation of Eq. ( TlUj) . One can rewrite 
Eq. ([2]) in the form 

N - ad 

d u q 



where 



1 N 

A = - ^ x i x jQij ~ ie ( A - 2 ) 

with given in Eq. ([6]). Integrating over q gives 

^ _ i-l) N m - D/2) f 1 ^ ^ *(£li - 1) 



T ° = (4^p7* J dxi "- dxN ±^dT2 • ( A - 3 ) 



The Landau equations for the representation (1A.3I) are [6] 

A = 0, 
= o 

dxi 



(A.4) 



Since A is a homogeneous function of Xj, the first equation in ( 1A.4j) is automatically 
satisfied when the second is. Eq. ( 1A.4I) is equivalent to Eq. ([5]), which means that the 



solution of Eq. ( 1A.4j) is an eigenvector of Q with zero eigenvalue. In general, Q has A" 
real eigenvalues Ai, . . . , A at. The characteristic equation of Q is given by 

/(A) = A^ + (-l^-iA^ 1 + (-1) V^A^ 2 - . . . (-l^VA + (-l) N a 

= (A-Ai)(A-A 2 )...(A-A n ) = 0. (A.5) 

For the case iV = 4 we have 

a = A1A2A3A4 = det(Q4), 

a\ = A1A2A3 + A1A2A4 + A1A3A4 + A2A3A4, 

a 2 = A1A2 + A1A3 + A1A4 + A2A3 + A2A4 + A3A4, 

a 3 = A 1 + A2 + A 3 + A4 = Tr(g 4 ), (A.6) 

Consider the case where Q has only one very small eigenvalue Ajv "C 1, then to a very 
good approximation 

X N ~ — , ai = AiA 2 • • • Atv_i ^ 0. (A. 7) 

<2i 
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Let V = {x®, ■ ■ ■ , x%} be the eigenvector corresponding to the eigenvalue Ajv- V is 
normalised to 

N 

5>? = 1. (A.8) 

t=i 

For latter use, we define 

v 2 = V.V. (A.9) 

The expansion of A around V reads 

1 N N 1 

A = - QijViVj + ^N^xlvi + 2 A ^ 2 _ ( A - 10 ) 

ij'=l i=l 

where yi = x\ — x® . In order to find the leading singularity, it will be sufficient to neglect 
the linear term in the rhs. The Q-matrix can be diagonalised by rotating the ?/-vector 

N 

y i = J2A ij z j , (A. II) 

i=i 

where A is an orthogonal matrix whose columns are the normalised eigenvectors of Q. 
Thus we have 

N spN I 

det ( A) = i, yj^ = ^ = -, 

j=l 

1 N ~ l 1 

^ = ^J2\Zi+^ N v 2 -te. (A.12) 

Note that the term \nz% in the rhs has been neglected as this term would give a con- 
tribution of the order 0(X 2 N ) to the final result. Eq. flA.3j) can now be re-written in the 
form 



N _ { -irT(N-D/2) r*~ 6(Y.*j ,,1,,,:,, 

° " 7_ ^■"^(Ef-A^ + A^-z.)^- (A ' 13) 



Although the original integration contour is some segment around the singular point Zi = 
with % = 1, . . . , N, the singular part will not be changed if we extend the integration 
contour to infinity, provided the power (N — D/2) of the denominator in Eq. ( 1A.13I) is 
sufficiently large. Integrating over z N gives 

T °" = ( \»I*2W D » 2)V L dZl '-' dZN ~ l ( E f = - Xi z 2 + W - ie)»-m ' (A ' 14) 

where the factor v comes from the 5-fuction. Asumming that Aj > for % — 1, . . . , K and 
Xj < for j = K + 1, . . . , N — 1 with < K < N — 1, we change the integration variables 
as follows 

U = \fkzi for i = 1, . . . K, 

t j = y/^TjZj for j = K + 1,...N-1. 1 ' ' 
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This makes sure that all are real. We get 
rpN _ (-l) N T(N-D/2)v 



n D/2 2 3D/2-N ^(-^N-K-l^ 

/ + OO /» + 00 
dti---dt K / dt K+1 ■ ■ ■ dtjv-i srn , (A. 16) 
7-oc + (-Efii^ + fc2 ) M/2 

where 

if 

b 2 = Y,^ + X nv 2 - ie, Re{b 2 ) > 0. (A. 17) 

i=i 

Changing to spherical coordinates and using the following formulae for the volume 

dt x ■■■dt K = / r K - x drdn K _ u / <m K -i = TV ^ /oV (A. 18) 

Jo J U A /2) 

we arrive at 

iN (-1) N T(N - D/2)v 2 7r( iV - A '- 1 )/ 2 



rpl 



n D/2 2 3D/2-N^^N-K-l ai T((N — K — l)/2) 

/+oo poo ^N-K-2 

oo dtl "' dtK J ^ \ b 2_ r 2 )N ~D/2 - ( A ' 19 ) 



Note that (b 2 - r 2 ) N ~ D / 2 = e -^N~D/2)^2 _ b 2^N~D/2 due tQ the fact that £>Q Uging 

[z+sy m 



"toj^-zMW-*™, (A . 20) 



gives 



T<f = — 1 lj 6 - — Tr^-^rffJV-D + g +l)/2) 

r D/2n3D/2-N . /? i W-it"-l. 



7T J/2 2 3D/2-JV ypTp. -,.-x ai 



dti---dt K = . (A.21) 

ii:f =1 t 2 + \ N v 2 -ie)( N - D + K ^)/ 2 

Repeat the above steps to write 

(^^(a^-D^ ^^((TV _ g + 1)/2) 
2 3 ^/ 2 "^ v /(-l) JV - A '- 1 a 1 (Atv^ 2 - ie)(^-D+i)/2 ' 1 ' ' 

This result was derived with the condition 

a x ^ and JV - £> + 1 > 0. (A.23) 

However it can be trivially analytically continued if we work in D = 4 — 2e so that it 
applies to iV < 3 in D = 4 by taking the limit e — > 0. 
Alternatively, with D = 4 and N = 3 the scalar function 

T = FT / ^1^2 77 a— T ^Tl 2 M ■ ( A ' 24 ) 

87r 2 y (Aizf + \ 2 z% + Aa^ - ie) 
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one first needs to dispose of the ultraviolet divergent first. To that effect we differentiate 
the above equation with respect to 77 = \ 3 v 2 with the result 



dT$ v 



dz\dz 2 - 



drj 8tt 2 7„ 00 (Aiz? + \ 2 z\ +V~ ie ) 2 

e in(2-K)/2 v 1 



(A.25) 



87r v /(-l) 2 -^AiA 2 V ~ ie' 
Integrating back (with respect to 77) we get 

in(2-K)/2 

^wpps'*" 2 -^^' (A ' 26) 

where C is a constant independent of 77. This result coincides with Eq. ( 1121 ). 



B Scalar box integrals with complex masses 

The derivation of the analytical expression of the scalar one-loop function for the box 
(N = 4) with complex internal masses in the most general case with no restriction on the 
external invariants is not tractable. However, if at least one of the invariant masses of 
the external legs is light-like one can derive an analytical formula in closed form starting 
from the the standard approach of 't Hooft and Veltman [28] (see also [37]). For our 
application there are at least 2 lightlike external momenta in all boxes. We explain here 
our derivation based on the method given in [2Sj for this special case. 

The scalar box integral is deduced from Eq. ( 1A.3I) with x 4 integrated out with the 
result 

D = (47r) 2 T 4 

1 

(ax 2 + by 2 + gz 2 + cxy + hxz + jyz + dx + ey + kz + f) 2 ' 

(B.27) 

where we have changed the integration variables as t = Ylt=i x ^ x = 12i=i x «> V = x i+ x 2, 
z = xi] and 

a = ^(<?33 + Q44 - 2Q 34 ) = pi, b = -(Q22 + Q33 - 2Q23) =pI, 

9 = ^(Qn + Q22 - 2Q 12 ) = pi, c = Q23 + Q34 - Q33 - Q24 = 2p 2 .p 3 , 

h = Q13 + Q24 ~ Qu ~ Q23 = 2pi.p 3 , j = Q12 + Q23 ~ Q22 ~ Q13 = 2pi.p 2 , 
d = Q34- Q44 = ml - ml - pi, e = Q 2A - Q 34 = m 2 2 - ml - p 2 2 - 2p 2 .p 3 , 

k = Qu - Q24 = m\-m 2 2 +p 2 l + 2pi.p4, f = - ie = ml - ie, (B.28) 

with Qij is defined in Eq. (ED- Our application will be to complex masses, m 2 , with 
? = 1,2, 3, 4, d, e, k, f are therefore complex parameters while other parameters are real. 
The two light-like external momenta can be either adjacent or opposite to each other. We 
consider in each of these tow cases separately. 



dx / dy 
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B.l Integral with two opposite lightlike external momenta 





Figure 17: A box diagram with two opposite lightlike external momenta p\ andpz- Double 
line means massless. 



For the box shown in Fig. [T7| with p\ = p\ = one gets a = g = and writes 



(13) 




dx I dy i dz 



Jo i^y 2 + CX U + h xz + 3U Z + dx + ey + kz + f) 2 
Integrating over z to get 



^B.29) 



D 



(13) 



y 



1 rx 

' dX J dy (Ax + B)(Cx + DY 



with 



A = cy + d, B = by 2 + ey + f 1 

C = (c + h)y + d, D = (b + j)y 2 + (e + k)y + f. 



One changes the integration order as 



1 px I>1 I>1 

dx dy = dy dx. 

JO JO Jy 



We get 



D o 13) = / d w I dx 

o Jy 



(Ax + B)(Cx + D) ' 



(B.30) 



(B.31) 



(B.32) 



(B.33) 
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where A, B, C, D are complex. Integrating over x as follows 

dx 



i 

dx 



(Ax + B)(Cx + D) AC J y ( x + f)(x+§) 



AD 




BC 




1 




AD 




BC 




1 




AD 




BC 




(B.34) 

where we have made sure that the arguments of the logarithms never cross the cut along 
the negative real axis. One easily gets 

D 0*> = f 1 d JL L A±± _ ln £±2\ m 35) 

J {cj - bh)y 2 + (dj + ck- eh)y + dk - fh \ Ay + B Cy + D J ! 

The discriminant of the quadratic function in the denominator of the prefactor is nothing 
but the Landau determinant. Indeed, 

det Q 4 = {dj + ck- eh) 2 - A(cj - bh)(dk - fh). (B.36) 

We write 

n (i3) _ 1 f 1 (I f, A + B _ C + D \ 

V ° ~ (cj-bh)(y 2 - yi )J [y-y 2 y-yJ^Ay + B Cy + DJ ' ^ ^ 

with 



-(dj + ck- eh) =f Vdet Q 4 

Vi = ; 7T\ • (B.38) 

y 2 2(cj -bh) y ' 

Now we have to look at the imaginary parts of the arguments of the logarithms in 
( IB. 371) . We write them explicitly 

A + B = by 2 + {c + e)y + d + f, 

Ay + B = (b + c)y 2 + (e + d)y + f, 

C + D = (b + j)y 2 + (e + k + c + h)y + d + f, 

Cy + D = (b + j + c + h)y 2 + (e + k + d)y + f. (B.39) 

Imaginary parts read 

Im(^4 + B) = lm(ey + d + /) = \m\ym 2 , + (1 — y)m\ — is] < 0, 

lvt\(Ay + B) = ha(ey + dy + f) — \m\ym 2 , + (1 — y)m 2 — ie] < 0, 

Im(C + D) = Im[(e + k)y + d + /] = lm[yml + (1 - y)m 2 3 - ie] < 0, 

lm(Cy + D) = Im[(e + k)y + dy + f] = lm[ym\ + (1 - y)m\ - ie] < 0. (B.40) 
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Using formula ln(a/6) = In a — In b for Im(a) Im(6) > 0, we rewrite f]B.37|) as 



r>(i3) 



with 



v/det 



Ax 

^3 

A, 



=F7f EEC" 1 )** f 1 dy-l—HA& 2 + B j y + C i ), (B.41) 



6 + c, B 1 = e + d, Ci = f, 

b, B 2 = c + e, C 2 = d + f, 

b + j, B 3 = e + k + c + h, C 3 = d + f, 

b + j + c + h, B 4 = e + k + d, C 4 = /. 



(B.42) 



We would like to make an important remark here. From Eq. ( IB. 401) we can re-write Eq. 
( IB. 371) in the form 



D 



(13) 



. A + B ^ Ay + B 
In — m 



(cj - bh)(y 2 - yi) J \y - y-i y-yij V C + D " Cy + D 
We notice that if y — y\ ]2 then AD = BC which means 



(B.43) 



A + B 



C + D 



Ay + B 

y=yi,2 Cy + D 



y=yi,2 



B 
~D 



y=yi,2 



[BAA) 



Thus, we get 



1 



o \y-V2 y-yi 



In 



A + B 



C + D 



Ay + B 



- In 

y=yi,2 Cy + D 



V=V1,2 



0. 



(B.45) 



Subtracting this zero contribution from Eq. ( IB. 411) we get another form 

2 4 



D, 



(13) 




EB-i: 



i+j 



^det(Q 4 ) i=1 j=1 



* dy HAjy 2 + Bjy + Cj) - \n(A 3 yf + B m + C 3 

'o y-yi 



which is more convenient for the evaluation in terms of Spence functions. 

Each integral in Eq. (1B.46|) can be written in terms of 4 Spence functions as given in 
Appendix B of [2H]- Thus L>J 13) can be written in terms of 32 Spence functions. 



B.2 Integral with two adjacent lightlike external momenta 

For the box shown in Fig. [18] with p\ = p\ = one gets b = g = and writes 



D 



(12) 



1 d r d r d 1 

Jo Jo ( ax2 + CX V + hxz + jyz + dx + ey + kz + f) 2 



[BA7) 
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Figure 18: A box diagram with two adjacent lightlike external momenta p\ andp 2 - Double 
line means massless. 



As in the case of D^ 3 \ integrating over z gives 



D 



(12) 



dx / dy 



dx / dy- 



-Ska^axy + b x 



(B.48) 



with 



s k = sign(Im(fc)), -s k a 1 = -s k (hx + jy + k) - ie', 
b\ = ax 2 + cxy + dx + ey + /, 
d\y + b\ = ax 2 + jy 2 + (c + h)xy + dx + (e + k)y + f — ie, 



(B.49) 



where we have used the fact that lm(aiy + b\) = lm[dx + (e + k)y + /] = Im[(x — y)m 2 + 
(1 — x)m 2 + ym\ — ie] < because < y < x < 1. 5 and e' are infinitesimal positive 
quantities which carry the sign of the imaginary parts of — Skdi and a x y + b\. For Jx, we 
integrate over y, similar to (1B.34[) . to get 



dy 



(ja — hc)y 2 + (jd — he — kc)y + jf — ke 



In 



(j + + k - ie' ^ (a + c)y 2 + (d + e)y + f 



hy + k — ie' 



ay 



dy + f 



Consider the prefactor 

det(Q 4 ) 

2/11(12) 



(jd -he- kc) 2 - A(ja - hc)(jf - ke), 
(he + kc- jd) =F ^det(Q 4 ) 
2(ja — he) 



[B.50) 



(B.51) 
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where the indices 11, 12 correspond to — and + signs respectively. We rewrite I\ as 
h 



1 Y(-v 



dy 

o y-vu 



' (j + h)y + fc - ie' _ (a + c)y 2 + {d + e)y + f 



hy + k — is' 



2 4 



- 1 VV(-l)^ / dy-^— 

V^tiQ^ttU Jo y-vu 



ay 2 + dy + f 



HA^ + B^y + C.j), (B.52) 



with 



A u = 0, B u = h, Cn = k, 

A 12 = 0, B12 — j + h, C12 = fc, 

= a + c, 5 13 = d + e, C13 = /, 

Au = a, B14 = d, C14 = /. 



(B.53) 



Thus l\ can be written in terms of 24 Spence functions. For I 2 we shift y = y + ax, a 
such that 



ja 2 + (c + + a = 0. 



(B.54) 



There are, in general, two values of a. The final result does not depend on which value 
of a we take. We have used this freedom to find bugs in the numerical calculation and 
it turns out to be a very powerful method to check the correctness of the imaginary part 
which can be very tricky for the case of equal masses. One gets 



"1 Ml—a)x 

I dx dy 

'0 J -ax 



1 



(Gx + H - ie')(Ex + F-is)' 



(B.55) 



with 



G = -s k h - s k ja, H = -s k jy - s k k, 

E = (2ja + c + h)y + d + a(e + k), F = jy 2 + (e + k)y + /. 



T3.56) 



For real a we have 



1 f(l—a)x 

dx / dy 

J —ax 



1 p(l—a)x rl p — ax 

dx dy — dx dy 

Jo Jo Jo 

1 — a /*1 p — a pi 

dy dx — dy dx. 

Jy/(l-a) Jo J-y/a 



(B.57) 



We write 



1 



1 



(Gx + H- ie')(Ex + F — is) GF - HE \Gx + H - is' Ex + F - is 



G 



E 



fB.58) 
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Integrating over x, we get 



ln G + H dy ln S + g 



a GF-HE E + F J GF -HE + F 



* ^ (B.5 9) 



The prefactor 
GF-HE 

Sk 

with 



o GF-HE Zy+F 



j(ja + c)y 2 + (2ajk + jd — he + kc)y + a{ke + k 2 — jf) + kd — hf 
j(ja + c)(y-y 21 )(y-y 22 ), (B.60) 



-(2ajk + jd-he + kc) + y/det(Q 4 ) , DC1 , 

2/21(22) = tpp — ' B - 61 

where the indices 21, 22 correspond to — and + signs respectively. We rewrite I 2 as 

h = — 1 Tj-iyiP, 

J- a y-y* E + F J y-y 2i ^ + F 

f~ a dy ^ + H 
+ / _^_i n ^ . (B.62) 

We make the substitutions y = y — a for the first integral, y = (1 — ot)y for the second 
integral and y = —ay for the third integral to get 



f 1 dy s k jy ~ s k h - s k k - ie' 

o y - a - y 2i jy 2 + (c + h + e + k)y + a + d + f-ie 
1 (1 - a)dy -s k (j + h)y - s k k - ie' 



o (l-a)y-y 2i (a + c + j + h)y 2 + (d + e + k)y + f - ie 

1 -ady -s k hy - s k k - ie 1 
In — — . (B.63) 

o -ocy - y2i ay 2 + dy + f - %e 

Consider the arguments of the three logarithms, as demonstrated in (1B.40I) . it is easy to 
see that the sign of the imaginary parts of the denominators is negative as indicated by 
— ie. The derivation is for real a. However, this result is also correct if a is complex as 
proven in [2H] • We can now rewrite I 2 as 

I 2 = 7 =L= £ f dy | - MA^y 2 + B 2j y + C 2j ), (B.64) 

s k Vdet(Q 4 ) Jo a jV - b j ~ Vn 
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with 



Ci 


= 1, ai = 1, b-i = a, 

7 i 7 ± 7 






Co 


= —(1 — a), a2 = 1 — a, 


62 = 0, 




Cq 
1 '& 


= — a, 03 = — a, 63 = 0, 






1 


= — 1 (I A — 1 ft /I = (T 






1 o 


= 1 — G!, CI5 = 1 — Oi, 65 


= 






= 0' (7,k = — O" ft« = 








= fioi = — sui Coi = 

w 5 -^zl ft: J 5 ^- zl 


— Suk — 


Si-h 


Y A 22 


w ; -^22 J K\J 1 ' 6 75 


n — 
^22 




^23 








^24 


= j, i?24 = c + /i + e + £;, 


C24 = 


a + d + f, 


^25 


= a + c + j + h, £?25 = 


+ e + fc, 




^26 


= a, i?26 = <^ C*26 = /■ 







(B.65) 

^2 can be written in terms of 36 Spence functions. Thus 

£>J 12) = l! + Sfc J 2 (B.66) 

contains 60 Spence functions. For the evaluation of in terms of Spence functions, it 
is better to do the following replacement for each logarithm in J 12 : 



\n(A 1:j y 2 + B ljV + Cy) -> \n(A ljy 2 + B ljV + Cy) - HA ljy 2 u + B ljVli + C : 
ln(A 2 y + B 2i y + C 2i ) -> ln(A 2j ?/ 2 + 5 2j -y + C 2i ) - ln(A 2i ^ + #2^2* + C; 



2j), 

(B.67) 



with i/2i = (y2i + bj)/dj. The argument for this is similar to that explained in the previous 
section, see Eq. ( 1B.46I) . 

For the boxes with one lightlike external momentum, the result is written in terms of 
72 Spence functions by using exactly the same method. 



C Singularities of the three point function 

In the main text we concentrated on the properties of the 4-point one-loop function es- 
pecially as concerns the occurrence of the leading Landau singularity which in that case 
is not integrable. Although a leading singularity in the 3-point function is integrable, it 
is instructive to study the case of the 3-point function in some detail as it sheds light on 
some properties we unravelled in the 4-point function. Moreover the three-point function 
appears also when shrinking or collapsing one of the internal lines into a point and there- 
fore its singularities are part of the singularities of the corresponding 4-point function. 
The study of the 3-point scalar integral is easier to handle as it involves less parameters. 
We take as an example, the 3-point loop integral shown in Fig. [19] that is part of the 
diagrams contributing to class (c). 
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Figure 19: A triangle diagram contributing to gg — > that can develop a leading Landau 
singularity for Mh > 2M\y and y/s~2 > m t + Mw, i-e. all the three particles in the loop 
can be simultaneously on-shell. 



In terms of the Passarino-Veltman appellation, this scalar integral writes 

T 3 (s 2 ) = C (s 2l Ml, 0, m 2 , M$y, M w ), s 2 = (p 4 + p 5 ) 2 . (C.68) 

The bottom-quark mass has been neglected by setting it to 0. The phase-space constraint 
on s 2 is M H < s 2 < s, see Eq. ( {TBI ). 

We will define the Landau determinant, detQs, corresponding to a 3-point function, 
Co(p 2 ,p 2 ,P3, m \i m 2, m i) according to the Passarino-Veltman notation 

detQ 3 (p 2 1 ,p 2 2 ,pl;m 2 1 ,m 2 2 ,ml) (C.69) 

with ?Tij the internal masses and p 2 the invariants of the external momenta. In the same 
spirit the determinant of a 2-point function will be defined as 

det Q 2 (p 2 ; ml, ml) = -\(p 2 ; mj, m 2 ,) = - ^p 2 - (m ± + m 2 ) 2 ^ ^p 2 - (m 1 - m 2 ) 2 ^(C70) 

where A(a, b, c) is the usual kinematic function, see Eq. ( fl5l) . For completeness and later 
reference, the determinant of the 1-point function is defined as 

detQi(m 2 ) = 2m 2 . (C.71) 

A necessary condition for a three point function to have a LLS is that it has exactly two 
cuts which can produce physical on-shell particles. The diagram in Fig. [Tj5] satisfies this 
condition when 



M H > 2M W and y/s^ > m t + M w . (C.72) 

These conditions are part of the conditions for our 4-point function (that we studied in 
section [3]) to have an LLS. In fact this three-point function is a reduced diagram from 
the point of view of our 4-point function where it is considered as a subleading Landau 
singularity. These conditions Eq. (1C.72I) represent the opening up of normal thresholds. 
We will refer to the first threshold M H > 2M W as the Higgs threshold (H — > W + W~), 
while the second condition will be referred to as the s 2 threshold (Hb — > Wt). 
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The sign condition (xj > 0), Eq. (jSj) for the case at hand are particularly simple here. For 
example, 

detg 13 = -M 2 H {m 2 + M 2 v )+2s 2 M 2 v <$, 

detQ 23 = -M 2 H (m 2 + M^) + s 2 (M 2 H - 2M^) < 0, (C.73) 
which together with Eq. ( 1C.72I) give 



a ^w^i M '' H - 2{m;+M ^ (a74) 



These inequalities are supplemented by the condition of vanishing Landau determinant 
in order for the appearance of the LLS. The Landau determinant in our case is 

det Q 3 (s 2 ,M 2 H ) = det(s 2 ,M£,0;m*,M^,A^) 



-2M^s 2 2 + 2M 2 H (m 2 + M^)s 2 - 2M 2 H M 2 H m 2 + (m 2 - M\ 



(C.75) 

We have chosen to pick up s 2 as the variable in which to study the location of the LLS, 
hence our notation det Qs(s 2 , Mfj). It is very rewarding to express this determinant in 
terms of a perfect square in s 2 plus a remainder which is the discriminant of the quadratic 
equation. We can then write 

detftt^) = -de tQl ( m? )(( 32 - 3 ^- J^ ?) ^g; ) with 

det Q 2 h = det Q 2 {M 2 H -M^M^) det Q 2 = det Q 2 (0; m 2 , M^) (C.76) 
4 = 2(m? + M^) + (M^-4M^)(l + ^i^). (C.77) 

It is important to note that the discriminant is the product of two sub-determiants, inde- 
pendent of s 2 , corresponding to two two-point functions each one obtained by collapsing 
or shrinking one of the internal lines bringing one vertex of the original 3-point function to 
coincide with the u s 2 vertex" , s 2 in which we write the perfect square. This is a general 
theorem [22] that applies to symmetric matrices based on the Jacobi ratio theorem for 
determinants [23] . 

The roots s 2t ± (from det Qs{s 2 , Mfj) = 0) give the position of the LLS as a function of 
Mjj, for fixed m t , Myy. In view of the constraint Eq. (1(1740 only one solution is possible. 
It is given by 

4 = 4 LS = (m 2 h {M 2 w + m 2 ) - (m 2 - M 2 W )M H ^/M 2 H - AM 2 ^ 

m 2 -M_ 2 

f 2 
Mr 



4 - \ M2 W M 2 HX /l - AM 2 W /M 2 H . (C.78) 



The surface that defines Eq. (10.780 is the surface of the LLS region. This surface is 
bounded however due to the constraint from the inequalities due to the normal thresholds 
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and the sign condition. This is what defines the region of the LLS singularity. In fact 
the normal thresholds are directly related to the range of the LLS region. First of all, if 

M 2 

M H < 2M W there is no LLS. At exactly the Higgs threshold, M H = 2M W and det Q 2 H = 
0, the LLS according to Eq. (10.781) occurs at S2 LS = 2[jn\ + Myy) which is the maximum 
value of S2 given by Eq. ( 1C.74I) . When Mh increases, the value of S2 LS decreases until S2 
reaches the S2 threshold, (m t + Mw) 2 , below which the LLS disappears. Therefore the 
S2 threshold, via the vanishing of the Landau determinant will give the maximum value 
of Mh for the appearance of the LLS. We therefore find that the region of the LLS is 
delimited as 

4M^ < M H < AM 2 W + M w )\ 

{m t + M w f < s 2 < 2{m 2 t +M v/ ). (C.79) 
Numerically, this corresponds to 

160.75GeV < M H < 172.89GeV, 

254.38GeV < < 271.06GeV. (C.80) 

This range in the variables Mh, S2 can be derived in a much simpler way. The Landau 
constraint of vanishing determinant detQs(s2, M H ) is a surface. This is bounded by 
tangents parallel to the coordinate variables [22], S2,M H so that with m t and My/ fixed, 
these extrema are given by 

d det Q 3 (s 2 ,M H ) = Q ^ g?t = s o ^ det Q f H = (gince det g 3(s2? M 2 } = Q) 
os 2 

=^ M H = 2M W =^ sf* 1 = 2(m\ + M^). (C.81) 
The other extrema are derived in a similar way by considering 

a det qg,, MS) = ^ AetQ? = ^ (c 82) 

s?"- 2 = (m, + M w f => M% = 4A4- + ^-(m, - M w f. 

m t 

It is crucially important to observe that these extrema do correspond to normal thresh- 
olds where a leading singularity and a sub-leading singularity coincide. This feature will 
be carried through to the case of the 4-point function. 

The location of the singularity, as well as its range, is well rendered in Fig. |2"U1 which 
shows how the location of the LLS moves as we vary the Higgs mass. Fig. [20] shows both 
the real and imaginary part of scalar 3-point function. Note that as shown in section [3], 
here the LLS is of a logarithm type. This explains why one observes a jump, a step 
function discontinuity, in the real part and a logarithmic singularity in the imaginary 
part or vice-versa. We see that for Mh = 159GeV < 2Myy, a funnel develops at the 
normal S2 threshold for the real part while the imaginary part develops a non-zero value 
past this threshold with a rather smooth and broad structure. For Mh = 2Mw, at the 
Higgs threshold, the imaginary part develops are very sharp dip at S2 = 2(m 2 +M v/ ) which 
is furthest from the normal S2 threshold at S2 = (m t + Mw) 2 - As the Higgs mass increases, 
this sharp dip moves to the left towards the normal S2 threshold beyond which the sharp 
peaks signalling the LLS disappear leaving only a dent at the normal (^2) threshold. 
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Figure 20: Left: the real part of Cq as a function of \fs~2 with various values of Mh- Right: 
the same plots for the imaginary part. 
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